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A model for the violation of the equivalence principle (VEP) on solar and reactor neutrinos is 
investigated. New limits for the VEP are obtained considering the mass-flavor mixing hypothesis 
and the VEP model. The search for an upper bound for VEP effects led to a possible solution where 
usual parameters tan 2 9 and Am 2 are practically unchanged. The set of best fit point parameters 
is given by \cj>Aj\ = 6.0±ll x 10~ 20 , tan 2 6 = 0.4781^ and Am 2 = 7.20±g;|f x 10" 5 eV 2 . This 
possible solution has increased the confidence level when compared with the MSW-only solution 
(tan 2 9 = 0, 462+£;o3g and Am 2 = 7, 75^q',i2 X 10~ 5 eV 2 ). A superior limit for VEP parameters has 
also been obtained: \$Aj\ < 3.9 x 10" 2O (3cr). 



I. INTRODUCTION 

Two decades ago, Gasperini introduced the idea of neu- 
trinos mixing between eigen-states of flavor and gravita- 
tional interaction, leading to a Violation of the Equiva- 
lence Principle (VEP) Q. The purpose of such model was 
to find a solution to the solar neutrino problem through 
an oscillation mechanism "a la" Pontecorvo supposing 
a minimal coupling between the gravitational field 4>{x) 
and a scalar approximation for the neutrino field. This 
approach just considers the neutrino kinetic energy con- 
tent as its mass when coupling with gravity. Later, Hal- 
prin and Leung 0, 0] introduced independently a cou- 
pled scalar approximated neutrino field with a linearized 
space-time metric, such that g^ v = rj^ v + where 
rjnv = diag(l, -1, -1, -f) and [J] = —2<p(x)S liV . 

Although these hypotheses have been formulated for 
massive neutrinos, no experimental data available at that 
time could distinguish between mass-flavor and gravity- 
flavor oscillations. It was much simpler therefore to con- 
sider only one of these two different effects. As a con- 
sequence, experimental confrontation made before the 
first KamLAND results Q considered this simple case 
of massless neutrinos, mixed only via gravitational inter- 
action (however see [g, 0, H| for a treatment with mass 
and VEP effects). In fact, a "just so" vacuum solution 
could explain all solar data. On the other hand, the 
increasing evidence of neutrino disappearance at short 
distances (= 180fcm) is not described by this kind of so- 
lution, which leads to neutrino oscillation lengths of the 
order of the Sun - Earth distance. 

With the increasing statistics on solar and reactor neu- 
trinos and the accumulated data from all other sources, 
one could ask what limits can be now imposed to neu- 
trino's VEP when we assume the mass-flavor mixing and 
MSW mechanism adding gravitational VEP interaction 
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in the system. In other words, would the neutrinos be 
good probes for VEP-like effects? (Other effects could 
behave with similar phenomenology and thus be indis- 
tinguishable from VEP, as will be discussed later). 

The VEP phenomenon manifests as a difference in the 
gravitational coupling for different states. In order to pa- 
rameterize its effects we will adopt the Post-Newtonian 
Parameterization [9(, where any difference from known 
gravitational Newtonian constant Gn is measured with 
a 7 factor, so that G' N = j m Gjy, where 7 m = j(m) de- 
pends on the mass to of the system. As the Gn constant 
is already been considered in the definition of the gravi- 
tational potential <p(r), one may also define the 7 factor 
as: 

4> = 7m , (!) 

where <fi is defined to be positive. For macroscopic bod- 
ies A and B, the difference between their ja and jb 
factors A7 = 7a — 7b has been measured with free 
fall experiments. Several gravitational sources are con- 
sidered as the Sun, Earth, and the galactic center, ob- 
taining A7 < 10~ 12 for astrophysical sources and 
A7 < I0~ 9 [ll| for terrestrial experiments. Interesting 
enough, some astrophysical events like pulsars with pecu- 
liar frequencies, could be explained if the neutrinos were 
experiencing VEP [12|, [l3|. By other hand, neutrinos 
cannot violate the equivalence principle by more than I 
part in 100 (90%C.L.) [l4| since one would be observing 
more erratic pulsars. 

II. VEP MODEL FOR MASSIVE NEUTRINOS 

We start stating that the model will apply only to 
weak gravitational fields, so that no spin-gravity effects 
will be considered here. By doing so, one may use the 
Klein-Gordon equation to describe the neutrino field: 
(<7 At! y9 M 9 I/ + to 2 ) ^ = 0, where g M „ is the metric tensor 
and \& represents a scalar approximation to the neutrino 
field. 

Following Halprin's approach, the metric tensor for a 
weak field can be written as g^ v = n^ + h fJ-u (x) where 
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ti mv = diag(l, -1, -1, -1) and = -2j m (f>(x)6 fiv [J], 
where the redefinition ([T]) is being used from now on. 
So, the Klein-Gordon equation with weak gravitational 
field is [(j7 M „ - 2 7m </)(x)(5 A1I ,) d^d" + m 2 ] * = 0. Assum- 
ing a plane-wave solution of the form ^ = ^oe l ^' x ^ Et \ 
one arrives at the energy-momentum relation for this in- 
teracting system: E 2 (l — 2 7m <^>) = p 2 (l + 2j m (f>) + to 2 . 
Using the fact that for neutrinos to <C p and ignoring 
terms with order higher than O(0 2 ), we finally have the 
energy-momentum relation for small masses and weak 
gravitational field: 



E = p (1 + 2 lm (j)) + — (1 + 4 7 „ 
2p 



(2) 



The above expression can be re-written as E = 

2 

E m + Eg so that E m = p + ^L. [ s th e free-particle 
energy-momentum relation (with to <C p) and E g = 
27m {p + is the gravitational contribution to the 
total energy. 

In order to introduce the neutrino mixing, one has to 
define basis on which each phenomenon takes place. The 
most general scheme for this model would be a three ba- 
sis system: a physical basis (states with definite mass), 
a weak base (states with definite flavor) and a gravita- 
tional base (states with definite gravitational couplings). 
This would mean that the dynamical and gravitational 
contributions to the total energy, E m and E g could not 
be simply added any more. Instead, the two physical 
quantities should be assigned to operators on different 
bases. Considering the further inclusion of weak interac- 
tions, and one third basis for it, the model will end with 
five free parameters Q (considering only two neutrino 
flavors). Although it is possible to carry on such analy- 
sis, it is interesting to test simpler models and, if any sign 
of VEP is found, a more complete analysis could be made 
in future works. To obtain a simpler model, we follow the 
hypothesis that the gravitational interaction takes place 
on the physical mass basis. This is exactly what has been 
done until now, when deriving the expression p]). 

Considering only two neutrino flavors, each eigen-state 
of mass has total energies E\ and E%, given by expres- 
sion ([2]), using m — * mi and j m — + 71 for E\, so that: 



Ex = p (l + 2 7l 0) + -i(l + 4 7l ^ 
2p 

and m —* mi and "/ m — > 72 for Ei\ 

2 

E 2 = p (1 + 2 72 0) + 2^.(1 + 4 72 ^ 
2p 



(3) 



(4) 



To describe a two level system, we introduce the Hamil- 
tonian 



AE f -1 Q 
2 V 1 



where AE = E% — E\ such that 



(5) 



AE = h 2p 0A 7 + - (7A771 2 + to 2 A 7 ) , (6) 

2p p 

where Am 2 = to 2 . — to 2 , A7 = 72—71, 7 = (72 + 7i)/2 
and fh 2 = (to 2 + m 2 )/2. 

Not all of these terms will contribute. Of the three 
terms with dependence on 1/p in ([5]), the last two (be- 
tween parentheses) are negligible, mainly because of the 
potential <j). Comparing all the sources of gravity that 
might have some effect here, as the Earth, the Sun, and 
larger scale structures such the Great Attractor 0, , 
the last contributes most, imposing a practically constant 
potential 4> ~ 3 x 10~ 5 [l7|, that is at least one order of 
magnitude more intense than the other sources Us- 
ing the definition of 7 in ([1]) and other VEP tests already 
cited, then 7=1 with A7 < 10~ 9 . These statements 
assure that (tAto 2 + m 2 A7) << Ato 2 , so that AE may 
be considered only as 



AE 



Am 2 
~2p~ 

Ag 
2E 



2p 4>A 1: 



(7) 
(8) 



where the usual consideration for neutrinos p = E was 
used and Aq = Ato 2 + 4E 2 (f>A'j is defined as an effective 
mass scale. 

We assume to-2 > mi. Nevertheless, the same relation 
does not have to hold for 71 and 72. Previous models 
for VEP considered only gravitational states for massless 
neutrinos. Consequently 7's could arbitrarily follow the 
hierarchy 72 > 71 in the same way it was done for the 
masses. In a model with VEP and massive neutrinos, the 
7's are dependent on the masses and so it is clear from the 
definition of Aq that the hierarchy between 71 and 72 will 
have influence on the resulting phenomenology. So we 
must consider two possibilities: if 72 > 71, following the 
same relationship defined for the masses, then A7 > 0; 
if 72 < 71 , then an inversed hierarchy on the VEP sector 
appears, and then A7 < 0. We consider 



A G = Ato 2 ± 4£ 2 



|0A 7 | 



(9) 



where |0A 7 | is one single parameter of the model and 
no further discussions about the single value of <f> are 
needed, as long as it is considered as a constant. The 
two possible relations between the 7's will be referred 
simply as +VEP and —VEP for the plus and minus 
sign on ©, respectively. 

When considering AE as a function of E, expression 
([7]) may vanish in the —VEP case. If so, some region of 
the oscillated spectrum might behave as the mass states 
were degenerate, even the mass eigen-states being differ- 
ent. For the +VEP case, AE can never be zero but it 
has a minimum value as its first derivative vanishes. Set- 
ting j^AE = for +VEP and AE = for -VEP, one 
obtains 
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E,, 



Am 2 



(10) 



which is then a critical energy of the model, either for 
+VEP and -VEP cases. 

As the energy E is a constant of motion, any previous 
solutions of the standard neutrino mixing model can be 
changed to accommodate the VEP hypothesis only by do- 
ing the substitution Am 2 — > Aq. Furthermore, no men- 
tion about the weak basis mixing was needed until this 
point. In this simplified two-bases version of the VEP 
model, gravity has no influence over the vacuum mixing 
(which would not be the case if a three-bases model is 
considered Q). 

The evolution of flavor states is given by the 
Schrodinger equation i4*^ = * (/) , with 



#(/) _ JJ $r(m) and 



(11) 
(12) 



where and represent the states and the Hamil- 
tonian, written on the flavor base (/). In general, states 
and operators in both bases are related through expres- 
sions pip and (|12p respectively, where U is a SU(2) 
transformation. When dealing with only two bases, U 
has only one physical parameter [H, HH which can be ex- 
pressed as 



U = 



cos 9 sin 9 
— sin 9 cos 9 



(13) 



In the flavor basis, one can introduce the effective weak 
potential 0[li]: 



V 



(/) _ ^ r _ 



W 



1 



(14) 



that describes the influence of a material medium on the 
neutrino conversion, known as the MSW effect [2(| HH • 
In expression Q14p. Gf is the Fermi constant and n e = 
n e (x) is the electron number density, used to describe 
globally neutral matter distributions. In our case, n e 
describes the Sun's and Earth's density profiles. The 

complete Hamiltonian is then given by = HW + 

( f) ~ 
V^r , where the sign denotes the presence of a mate- 
rial medium. Plane wave solutions of the Schrodinger 
equation for this system are of the form ^^\t) = 

exp "to . Considering the initial state to 

be constituted only by electron ( anti- ) neutrinos [38j . one 
arrives at the survival probability P ee (t) = |\&g\I/(-^(t)| 2 , 
which describes the probability that a neutrino is both 
produced and detected as an electron (anti-) neutrino. 
The simplest solution corresponds to the vacuum case 
(VAC), where n e = 0, and is given by: 



1], (L.El = 1- si.r 20 siir ( 



(15) 
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(b) VAC- VEP case. 
FIG. 1: Survival probability P ee (E) for VACiVEP. Figure 

Hal corresponds to VAC+VEP andflblto VAC- VEP. The 
parameters corresponding to the usual neutrino mixing are taken 
to be sin 2 26 = 0, 86, Am 2 = 8, X 10~ 5 eV 2 and L = 180fcm ^ 
(this last one corresponds to the average distances considered for 
KamLAND). Each line represents a specific value of VEP: 
\<t> A 7 | = (no VEP), 1CT 21 (Ei, = UOMeV) and IO" 20 
(B* = 45MeV). 



where L is the distance between the source and the de- 
tector. The resulting periodic pattern has an oscillation 
length A = AitE /\Ag\, where the absolute value of Aq is 
used since it can become negative (in the —VEP case). 
In practice, L is fixed (a characteristic of the experiment) 
and we observe P ee as a function of E. If VEP is not 
present, A depends linearly on E and the oscillation stops 
when E S> LAm 2 /4tt. Otherwise, the dependence of Aq 
on E prevents the oscillation from stopping as now A has 
a maximum where \Aq\/E has a minimum, at E — E+. 
In a particular case, A — > oo when E — > E+, for the 
—VEP scenario. 

Figs. Hal and [TBI show expression (TT5|) for +VEP and 
—VEP respectively. The values used for sin 2 29 and Am 2 
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are those found in the literature [23| for the standard mix- 
ing, and L is constant and refers to the KamLAND [22j 
experiment average source-detector distance. On both 
figures, one can observe the new effect where the oscilla- 
tions are restored for energies above E+. The presented 
values of \(j> A7I are chosen so that any predicted new ef- 
fect will not be visible within the reactors' spectra energy 
range (approximately E < 9MeV) . This give us a visual 
"first limit" for VEP as \<f>A-y\ < 1CT 20 , if the present 
data reject the hypothesis. 

Solutions that describe solar neutrinos must consider 
the Sun's matter profile, given by the Solar Standard 
Model (SSM) The MSW effect predicts not only 

conversion between flavor states but also, under certain 
conditions, conversion of mass states referred to as non- 
adiabatic effects [l8l . [l9j . To better understand these ef- 
fects, one has to transform H^' back to the mass ba- 
sis, where it should be diagonal. The introduction of the 
weak potential V~w assures that this transformation is dif- 
ferent from U. It's then necessary to define the effective 
mass basis so that §( m > = l/t#W and = U^H^U 

where ff( m ) has a diagonal form. The transformation U 
is defined as (fT5)) with 9^9. Requiring ij( m ) to be 
diagonal, one arrives at the effective mixing in matter, 
given by 



cos2#(x) = 



A G cos 29 - A(x) 



[A G cos 26 - A{x)\ + A G sin 2 9 



. (16) 



where A(x) = ^/2GpEn e (x). The Schrodinger equa- 
tion will not retain the same form under such transfor- 
mation, since U has a dependence on the position x. 
Transforming states and the Hamiltonian from the fla- 
vor to the effective mass basis results in i U^-^U vj/M = 
if( m )\lH m ). The resulting evolution operator has addi- 
tional off-diagonal terms that come from the derivative 
U^-j^U, which together with the diagonal Hamiltonian 
fji m ) give us 
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(a) MSW+VEP case. 
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(b) MSW- VEP case. 

FIG. 2: Adiabaticity Coefficient F(E) for MSW+VEP. 
Figure l2al corresponds to MSW+VEP and|2E]to MSW- VEP. The 
parameters corresponding to the usual neutrino mixing are taken 
to be sin 2 26* = 0, 86, Am 2 = 8, X 10" 

profile is the one from BS05(OP) [I 

specific value of VEP: \<j> A-y| = (no VEP), 10" 20 , lO" 19 , 
10 — 18 and 10 -17 . For the +VEP case, the system is adiabatic 
For —VEP, non-adiabatic effects occur when E — > E+ for any 
value of |0 A7I 



5 eV 2 and the Sun's matter 
J. Each line represents a 
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where E\ and E2 are the eigen-values of 7j( m ) and 9 is 
implicitly given by (fT6|) . The off-diagonal terms in (fT7|) 
result in a non-zero probability of conversion between ef- 
fective mass states. The intensity of these non-adiabatic 
effects can be measured by the relation between the diag- 
onal and the off-diagonal terms of (fTT)) . such that when 



En, — Ei 



the evolution operator on (jTTJ) is ap- 
proximately diagonal. This condition can be summarized 
in the form of a Adiabaticity Coefficient T(x, E) defined 



dx 



AE 



< 1 , with AE = E 2 - Ex , (18) 



T(x,E) 



where AE = [A G cos 20 - A(x)f + A 2 G sin 2 29. 
To better appreciate non-adiabatic effects on the neu- 
trino spectrum, it is useful to define T as a function 
of the energy only, eliminating the x dependence by 
taking the maximum value of T(x,E) for any E, i.e. 
T{E) = max{T(x,E)}. Fig. [2a| shows T(E) for MSW 
only in comparison with the MSW+VEP case, for some 
values of |0A7|. Inside the solar neutrino spectrum re- 
gion T(E) is never higher than 10" 3 and so non-adiabatic 
effects are not expected in this case. On the other hand, 
fig. [2b] reveals that extremely non-adiabatic effects occur 
in the characteristic energies E+ for the MSW— VEP case. 
Such behavior happens when A G — > 0, causing both E\ 
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and E2 to vanish. As a consequence, the off-diagonal 
terms of (fTT)) become infinitely larger than the diagonal 

ones (as these goes to zero), even when 4^ is naturally 
small, as they are expected to be in the Sun (as can be 
seen from T(E) in figs. I2"al and [2b|l . 

For those cases where condition (|T8|) holds, equation 
(|17p may be solved in the adiabatic approximation that 
leads to the following survival probability [HI : 



1 + cos26» cos26»(a;) 



sin2# sin20(a;) cosa(a;) 



(19) 



where cos2# is also given by (|16|) where cos26> = 
cos20(xo) being xq the neutrino production point. The 
factor cosa(a;) corresponds to the oscillating term with 
a(x) — Jq AE(x') dx' . When the matter contribution 
vanishes, expression (|19j) corresponds to vacuum solution 
([15|) . Moreover, if Aq ^> 10~ 10 eV 2 , cosa(x) rapidly os- 
cillates (when compared to the Sun's dimensions [l^.[25l|) 
and is ruled out by the average over the production 
point xq. Without VEP this condition is satisfied as 
Am 2 = 8.0 x 10- 5 eV 2 [Hj], what leads to the useful sim- 
plified survival probability for solar neutrinos, 



P ee( x ) = 7:[l + cos26l cos26»(x) 



(20) 



From the studies of the adiabaticity coefficient, the 
above expression is expected to hold for the MSW + 
VEP case and almost everywhere for MSW — VEP, 
except in the neighborhood of E+. Fig. [3al shows a com- 
parison between expression (|20p with no VEP and with 
\<f>Ay\ = 5 x 10- 20 , for the MSW + VEP case. This 
value of VEP was chosen so that E+ is just above the 
solar neutrino spectrum (E+ = 20MeV). The survival 
probability with VEP is always greater than the one for 
MSW only, but this difference only becomes appreciable 
for energies above E*. Fig. I3bl shows the same compar- 
ison for MSW— VEP. As expected, non-adiabatic effects 
occurs near £7+. Fig. I3bl also shows a numerical solution 
of the equation (|17p , in which the non-adiabatic behavior 
can be seen in details. These effects are confined inside 
a narrow region around E+ and they are not observable 
with the present data statistics. The adiabatic approx- 
imation (j2"0|) describes well the survival probability by 
any practical means, in both ±VEP cases. 

When arriving at night time, solar neutrinos cross sev- 
eral Earth layers. Again, the presence of matter alters 
the survival probability in a way that night time neutri- 
nos have more chance to survive than those arriving at 
day time. This effect is called regeneration and is not 
observed on the solar neutrino data |26[. On the same 
way that new non-adiabatic effects were predicted for 
day neutrinos in the MSW — VEP case, new regenera- 
tion signal may also be expected. To account for theses 
possibilities a numerical solution of (TIT]) , for the Earth's 
matter profile, is used. Fig. 2] is the equivalent to [3bl 
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FIG. 3: Survival probability P ee {E) for MSWiVEP. Figure 

l3al corresponds to MSW+VEP and [3b] to MSW- VEP. The 
parameters corresponding to the usual neutrino mixing are taken 
to be sin 2 29 = 0, 86, Am 2 = 8, X lCr 5 eV 2 . On fig[3al 
expression H20H compares the case without VEP and for 
|r/>A7| = 5 X lCr 20 . Fig. l3alalso compares a numerical solution 
for 117> that shows new non-adiabatics effects in the proximities 
of Ei, (in detail). 



for night neutrinos. In the neighborhood of E+, the solar 
non-adiabatic effects are intensified by Earth's matter. 
Fig. [S] shows a measure of the asymmetry between night 
and day for the —VEP case. As it can be seen, an ex- 
cess is expected in a region wider than the one where the 
solar non-adiabatic effect takes place. The absence of re- 
generations signs on the solar neutrinos data imposes a 
stronger limit on E+ for the —VEP case than for +VEP. 
This has direct consequences on the limits for |0A7|. 



III. DATA ANALYSIS 

For solar neutrinos, we consider data from 
Homestake H3, Sage [H|, Gallex/GNO [H[, Su- 
perKamikande(SK) [H and SNO (I [H| e II (H) 
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where A*' 1 and Nf x are the theoretical and experimental 
countings. 

The x 2 distribution is taken to be a function of the 
three parameters: tan 2 8, Am 2 and | ^> A.'y] . In order to 
identify local minima, we minimize \ 2 with respect to 
two parameters, varing the third freely. This procedure is 
repeated for each parameter, for both +VEP and —VEP 
cases. Figure H] shows A\ 2 as a function of |</>A7|, where 



FIG. 4: Survival probability P ee (E) for neutrinos 
arriving AT NIGHT, WITH MSW— VEP. A comparison between 
day and night probabilities shows an excess for the night time. 
The parameters used in this plot are the same as in fig. I3bl 
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FIG. 5: Day-Night asymmetry for MSW-VEP. The -VEP 

hypothesis predicts new regeneration effects due to Earth's 
matter generating an excess of solar neutrinos arriving at night 
time, for energies close to E+. 
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experiments. New data from SNO-III phase [33| and 
Borexino [34| were not included, but we don't expect 
them to change our conclusions. As for reactor anti- 
neutrinos, KamLAND data is considered (22J. A x 2 
analysis is done, where we define 

X 2 =xlun+X 2 KL- (21) 

The solar neutrinos contributions xtun * s given by (26l |: 

119 

xLn = E [Rt°-#r p ] [s 2 ];-[Rf-RT] ■ ( 22 ) 

where Rf 1 and Rf° are the theoretical and experimental 
rates respectively and S 2 takes in account all the corre- 
lation between uncertainties [32l . |35| . Reactor neutrinos 
contribute through a Poisson statistics (36j : 



FIG. 6: Ax 2 DISTRIBUTION AS A FUNCTION OF </>A7|. The 
dashed curve represents AXs Un , while the dotted one refers to 
Ax|f r. The continuous curve stands for the global fit. Figures l6al 
and I6bl shows the +VEP and —VEP cases respectively. 



When |0A7| < 10~ 22 , the analysis shows that the 
standard global solution for solar and reactor neutri- 
nos is recovered: tan 2 6 = 0. 462+^35 and Am 2 — 

?-75±o.i2 x 10~ 5 eF 2 . 

For the +VEP case, figure l6al shows a global minimum 
with better fit than the no- VEP case. Both solar and 
reactor neutrino alone shows a minimum, although in the 
solar case this is less than one sigma. So it is better suited 
to say that the solar data excludes +VEP with I0A7I > 
3.4 x 10~ 19 (3<t). On the other hand, the fit with reactor 
neutrino data from KamLAND is enhanced by almost 
3(7 what is, to say the least, surprising. The VAC + 
VEP solution seems to almost exclude VAC alone by a 
3cr factor (A% 2 w 8). This difference is smothed when 
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L Q /E v (km/MeV) 



FIG. 7: Comparison between VAC + VEP, VAC-alone and 
KamLAND DATA. There is an overall reduction in the distance 
between the VAC + VEP fitted model and the data, in 
comparison with the VAC-alone fit. 



both solar and reactor neutrinos are taken in account 
together. In this case, the MSW + VEP solution fits 
the dataset better than MSW alone by a difference of 
Ax 2 ~ 5, with the best fit point parameters given by 



+3.1 
-2.3 



XlO 



-20 



tan 2 9 = 0.461_f g;^ and Am' 2 



|(/>A 7 | = 6.0 

7.20i~o;27 x W~ 5 eV 2 (errors corresponding to ler). These 
values are just la apart from those corresponding to the 
standard MSW fit. 

In order to appreciate and justify the improvement we 
obtained in the fit of KamLAND data with VAC + VEP 
solution when compared with the VAC alone fit, we in- 
troduce figure El Although it is generally accepted that 
VAC alone generates a good fit of these experimental 
data, one can see from this figure that the introduction 
of VEP effects makes an additional improvement of the 
quality of the fit, especially in the higher energy region of 
the reactor spectrum. The VAC-alone model predicts a 
larger "fall" in the left side than indicated by data, while 
the VAC+VEP solution seems to be in better agreement 
in this region. This can also be seen in figure [H which is 
similar to figure [Taj but now shows the VAC solution in 
comparison with the just found VAC+VEP for the b.f.p 
of global solution. 

As for the —VEP case, figure [6b] shows no local min- 
imum in the global Ax 2 - When only the solar data is 
taken in account, two minima arc visible although not 
so evident. They correspond to |</>A7| = 9.0 x 10~ 20 
and |</>A7| = 2.0 x 10~ 19 . The mixing angle is roughly 
the same but the mass scales are Am 2 = 6.2 x 10~ 5 eV 2 
and Am 2 = 1.6 X 10~ 4 eV 2 respectively. These values 
correspond to the same critical energy: E* = \AMeV 
(within the la range). But this is not consistent with 
the KamLAND data, which clearly excludes — VEP with 
|0A7| > 3.9 x 10~ 20 (3a). So this can be regarded 
as a superior limit, strongly imposed by KamLAND 
data. If one discards for a moment the best fit point 
of the +VEP case, a superior limit is also obtained: 




-+- 



4- 



-VAC 
- VAC+VEP 

H 1 r 



123456789 10 
Neutrino Energy - E (MeV) 

FIG. 8: Comparison between VAC + VEP and VAC alone 
SOLUTIONS. The parameters used in this plot are given on table [I] 
For the lower energies, the difference is small, but high energies 
neutrinos are expected to have a greater survival probability. 

Parameter MSW/VAC MSW+VEF 



|0A 7 j(xlO— ' ) 
Am 2 (xl0- 5 eV 2 ) 
tan 2 6 



7.75 
0.462 



zero 

0.16 

0.12 

0.1)43 
0.03<i 



6.0 
7.20 
0.478 



T3"T 

2.3 
+0.35 
0.27 
+0.040 
0.038 



TABLE I: Possible solution for the VEP model in solar and 
reactor neutrinos. The second row shows the MSW-only case, 
where the standard mass-flavor solution is recovered. The third 
shows the VEP solution fot MSW+VEP. 



I0A7I < 1.5 x 10 19 (3<t). The superior limit for pres- 
ence of VEP represented by the M SW ± VEP model 
(regardless of the sign case) is the less restrictive of these 
values, so that: 



\4>A-y\ < 1.5 x 10~ 19 (3cr), 



(25) 



since the + and —VEP cases are obviously mutually ex- 
cluding. 

In order to compare with the macroscopic experiments 
on VEP, one has to consider an estimative of the grav- 
itational potential </>. It seems that, among several 
possible sources, the Great Attractor offers the largest 
contribution Q, with its best estimative given by <j) — 
3 x 10~ 5 . So the upper bound of |<^>A 7 | given in (|2"5")l cor- 
responds to a maximum value of the order IA7I < 10~ , 
wich is in agreement with the macroscopic experiments 
presented in the introduction (IA7I < 10~ 12 ). The possi- 
ble solution found with non-zero VEP value corresponds 

tO |A 7 | RS 10~ 15 . 



IV. CONCLUSION 

The results of our analysis suggest that VEP-like ef- 
fects might be taking place in the neutrino sector. The 
model offers two theoretical possibilities: one in which 
greater mass represents greater gravitational coupling 
(here called VAC/M SW + VEP) and an inverse sit- 
uation, where greater mass implicates a smaller cou- 
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pling with the gravitational field (VAC/MSW - VEP). 
With the latest statistics presented by the KamLAND 
collaboration and all known solar neutrino data, ex- 
cept Borexino and SNO-III, such evidence manifests it- 
self as a surprising enhancement of the global fit for the 
VAC/MSW + VEP case. The analysis showed that the 
VAC model alone is almost excluded with 3a. In the case 
of reactor neutrinos, the VAC/MSW + VEP model pre- 
dicts a higher survival probability for high energy reactor 
neutrinos, which is shown to be in better agreement with 
the experimental data. The VAC/MSW - VEP option 
is excluded by KamLAND data within the present statis- 
tics. All results are summarized on table |TJ 

These evidences are enough to encourage further 
studying. Hopefully future results from the KamLAND 
collaboration, with increased statistics, would allow not 
only a better analysis of this version of the model but of 
more complex versions as well. In particular, this anal- 
ysis shows that models with phenomenology similar to 
this one might enrich the neutrino sector of the Stan- 
dard Model. The VEP hypothesis presented here is just 



one. Any model that presents a mixing scenario, with a 
Hamiltonian like ([5J and with AE given by an expression 
with the same momentum dependency as the one seen on 
([7]), would lead to a better fit. As an example, the com- 
bination of Violation of Lorentz Invariance (VLI) models 
with mass-flavor mixing presents the same phenomeno- 
logical behavior (3?| as shown on §5§ and ([7]), needing 
only a parameter reinterpretation: \a\ = 2|0A7|, where 
| a | 7^ implies VLI between neutrino flavors (as defined 
on section 2 of [37]). In general, the results of table U also 
apply to a possible Lorentz violation or a combination of 
both phenomena. 

Considering only VEP effects, this analysis has im- 
proved the previous limits for A7 in two orders of mag- 
netude, lowering it from IA7I < 10~ 12 down to | A.'y | < 
10~ 14 , when it is assumed that neutrinos are mainly af- 
fected by the gravity of the Great Attractor. 

We would like to thank FAPESP, CNPq and CAPES 
for several financial supports. 
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